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1. Introduction 

Renormalization of lattice operators is an essential component for many precision applica- 
tions of lattice QCD. Amongst these are the determination of quark masses and weak interaction 
phenomenology, both of which require matching operators to continuum renormalization schemes. 

Non-perturbative renormalization (NPR) [|l|] is a technique for performing such matching cal- 
culations using lattice simulations directly, and does not require the use of lattice perturbation the- 
ory. Instead, it requires that the renormalization scale jjL is well separated from both the cutoff scale 
of the simulation, \/a, and the scale at which non-perturbative QCD effects become significant: 

Aqcd<jU<1/«- (1-1) 

There are relatively few NPR studies using staggered fermions [Q, ^, ^, ^, ^. We have calcu- 
lated matching factors for all staggered fermion bilinears spread over a 2^ hypercube, determining 
their renormalization scale dependence and comparing them to the predictions of one-loop lattice 
perturbation theory (PT) [^]. We use both HYP-smeared and asqtad valence fermions on an asqtad 
sea. Previously we used the MILC coarse ensembles ||5|]; here, in Section ^ we extend the results 
to the fine lattices. 

The bilinear operators form the building blocks for the four-quark operators used for weak in- 
teraction phenomenology. An ongoing computation of Bk using HYP-smeared staggered fermions 
finds Bk = 0.727 ib0.004(stat) ib0.038(sys) The dominant error comes from using the one- 
loop perturbative matching factor, and is estimated to be 4.4% assuming that two-loop effects are 
of size a?. The present work is an important first step towards fully non-perturbative matching, 
which will both reduce the error and improve its reliability. The current study will give a general 
indication of how accurate one-loop predictions are, and thus indicate the reliability of the error 
estimate used in [j8|]. The bilinears are also of phenomenological interest in their own right for the 
determination of quark masses. 

The asqtad light-quark matching factor was calculated non-perturbatively in |j3|], using the 
RI/MOM scheme. It is however advantageous to use a matching scheme that is less sensitive to the 
effects of low-energy QCD. In Section || we present results for the scalar and pseudoscalar bilinear 
channels using the RI/SMOM scheme. We find a significant reduction in non-perturbative splitting 
between these channels and weaker mass dependence relative to the RI/MOM case. 

2. Exceptional NPR 

We compute bilinear Z-factors for all spin (/s) and taste (i§f) structures. These fall into 35 
irreducible representations (irreps) under the lattice symmetry group, which are enumerated in 
Table Unlike with Wilson or domain wall fermions, staggered bilinears are slightly non-local, 
spread over a 2^ hypercube and requiring gauge links in their definition. Their connection with the 
continuum bilinears may be expressed as follows: 

eW(75 ® ^f)Q'{x) ^ Z^^^^^ ^^^(n) {ys<^^F)ABUn+A,n+B M ■ (2.1) 

A,B 

Here Xa{}T') is a single-component field located at site n+A, where n labels the hypercube and 
A a site within the hypercube. The Q{x) are four-component Dirac spinors that also come in 
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Table 1: Covariant bilinears forming irreps of the lattice symmetry group. Indices /i, v and p are 
summed from 1—4, except that all are different. Pseudoscalar and axial bilinears are not listed: 
they can be obtained from scalar and vector, respectively, by multiplication by 75 (8) 1^5. Bilinears 
related in this way have the same matching factors. This operation also implies the identity of the 
Z-factors for the tensor bilinears within square brackets. Bilinears mai^ked in blue are used as the 
denominators of the ratios discussed in the text. 



four "tastes" (indices left implicit). Un+A,n+B is an average over all minimal length gauge paths 
connecting n+Aton + B. 

We compute quark propagators on coarse (a 0.12 fm) and fine (a w 0.09 fm) MILC gauge 
ensembles, using both asqtad and HYP-smeared valence actions at three masses for which am^ai = 
anisea- The use of momentum sources results in small statistical eiTors with relatively few (~ 10) 
configurations. One feature that distinguishes NPR using staggered fermions is that for each phys- 
ical momentum p, 16 inversions must be performed corresponding to p + nA and the resulting 
momentum-space propagator is a 16 x 16 matrix. 

The gauge links used in the the construction of the bilinears are smeared. For the HYP valence 
action we use HYP smearing; for the asqtad valence action we use "Fat 7 -1- Lepage" smearing. 

We impose RI7M0M renormalization conditions on the bilinears of Equation ( p] ) for a num- 
ber of physical momenta (~ 10) at each mass, and extrapolate these results to the chiral limit. The 
pseudoscalar channel has a non-perturbative l/m dependence and we do not present the results for 
it here. 

We present our comparisons to PT in terms of ratios of Z-factors with the same spin-structure 
but different tastes. This makes the perturbative predictions especially simple because the contin- 
uum running cancels in the ratio: 

Zsi»Fl jp) _ . «(Mo) r^LAT ^LAT 1 

- 1 + ^ [(^S^F2 - Csr^Fll , (2.2) 

where ^Uq is on the order of the lattice scale and is not related to the renormalization scale p. Since 
the ratios calculated non-perturbatively are sensitive to lattice artifacts and low-energy effects of 
QCD, their variation with respect to p provides an indicator as to whether one is in the NPR 
window. 

Results for the fine HYP and asqtad ratios at a single scale p = 2.1 GeV {{apY = 0.81) are 
presented in Figure |l| and Figure ^ The perturbative results aix evaluated with a (3 GeV) = 0.24. 
For the HYP ratios, we see good agreement with PT for the vector, axial vector, and tensor ratios at 
around the percent level. The scalar channel shows agreement at the 5 - 10 % level. Similar results 
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Figure 1 : Comparison of Z-factor ratios for vector, axial, tensor and scalar bilinears to perturbation 
theory for HYP fermions on fine MILC lattices. Horizontal lines show perturbative predictions, 
with solid/dotted lines showing results with/without mean-field improvement. Results are in the 
chiral limit for the momentum described in the text. 



were found on the coarse MILC lattices [^. Overall, they indicate agreement at or better than the 
naive expectation of &{a?) truncation errors. The variation of these quantities with renormalization 
scale is shown in the left-hand column of Figure ^. 

For the asqtad ratios, the vector, axial vector, and tensor agree with perturbation theory at the 
few-percent level, while the scalar ratios show significant discrepancy (~ 25%). The variation of 
these quantities with renormalization scale is shown in the right-hand column of Figure ||. For both 
the HYP and asqtad cases, the nonperturbative values in the scalar channel more nearly approach 
the perturbative predictions as (apf- increases. 



3. Non-exceptional NPR 

Here we present our first results for staggered NPR using a non-exceptional scheme (RI/S- 
MOM) [Q]. Non-exceptional schemes have reduced sensitivity to nonperturbative effects com- 
pared to their exceptional counterparts, increasing the range of the NPR window. They achieve this 
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Figure 2: As for Fig. |l| except for asqtad fermions. Only results from mean-field improved PT are 
shown. 



by using a kinematic setup that retains a single kinematic invariant but for which no channel has 
zero momentum flow (i.e. there are no "exceptional" channels). In the case of bilinears, the flow 
through the external quark lines satisfies p\= p\ = {p\ — PiY- Figure ^ shows the nonperturba- 
tive splitting between the scalar and pseudoscalar channels in the exceptional and non-exceptional 
cases. Figure ^ compares the mass dependence of the scalar channel for the exceptional and 
non-exceptional cases. The non-exceptional setup has a much weaker mass dependence. This will 
improve the determination of using the asqtad action, where the strange sea-quark mass is fixed 
and cannot be extrapolated to zero. 



4. Conclusion 

We have presented non-perturbative matching factors for staggered bilinears using both HYP- 
smeared and asqtad valence actions, comparing our results to the predictions of one-loop pertur- 
bation theory over a range of momenta. We generally find good agreement (up to ^(a^)) for the 
Z-ratios presented here, except in the asqtad scalar- channel (particularly for small {apY). 

We also present the first study of non-exceptional schemes using staggered fermions. We ver- 
ify that nonperturbative splitting between the scalar and pseudoscalar channels is greatly reduced 
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HYP scalar ratios (fine) 



ASQ scalar ratios (fine) 



1.45 
1.40 
1.35 
1.30 
1.25 
1.20 
1.15 
1.10 
1.05 
1.00 



0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

[apY 



1.55 
1.50 
1.45 
1.40 
1.35 
1.30 
1.25 
1.20 
1.15 



0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

{apf 



Figure 3: Momentum dependence of Z-ratios for HYP (left column) and asqtad (right column) 
fermions, calculated on fine MILC lattices. The colors correspond with those of the data plotted at 
a single scale in Figure |I] (HYP) and Figure ^ (asqtad). 
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(a) Non-perturbative splitting of the scalar and pseu- (b) Mass dependence of the scalar channel in the 
doscalar channels using an exceptional (E) and exceptional (E) and non-exceptional (NE) schemes, 
non-exceptional (NE) scheme. 



relative to the exceptional case. The mass dependence of the scalar channel is also much milder, 
which should lead to an improvement in the determination of using this action. 

Detailed descriptions of these results and associated theoretical work are in preparation [|T0|]. 
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